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Abstract 



New families of exact general relativistic thick disks are constructed 
using the "displace, cut, fill and reflect" method. A class of functions 
used to "fill" the disks is derived imposing conditions on the first and 
second derivatives to generate physically acceptable disks. The anal- 
^jr^' ysis of the function's curvature further restrict the ranges of the free 

parameters that allow phisically acceptable disks. Then this class of 
functions together with the Schwarzschild metric is employed to con- 
. struct thick disks in isotropic, Weyl and Schwarzschild canonical co- 

' ordinates. In these last coordinates an additional function must be 

added to one of the metric coefficients to generate exact disks. Disks 
in isotropic and Weyl coordinates satisfy all energy conditions, but 
those in Schwarzschild canonical coordinates do not satisfy the domi- 
nant energy condition. 

PACS numbers: 04.20.Jb, 04.40.-b 



1 Introduction 

Exact solutions of Einstein's field equations with axial symmetry play an 
important role in the astrophysical applications of general relativity, since 
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the natural shape of an isolated self-gravitating fluid is axially symmetric. 
In particular, disk-like configurations of matter are of great interest, since 
they can be used as models of galaxies and accretion disks. 

Solutions for static thin disks without radial pressure were first studied 
by Bonnor and Sackfield [Q, and Morgan and Morgan and with radial 
pressure by Morgan and Morgan "T. Several classes of exact solutions of 
the Einstein field equations corresponding to static thin disks with or with- 
out radial pressure have been obtained by different authors |l]-jl3j. Thin 
rotating disks that can be considered as a source of the Kerr metric were 
presented in |^, while rotating disks with heat flow were studied in [T^ . 
Also thin disks with radial tension JH]) magnetic fields |E1 and magnetic 
and electric fields jTQ were considered. The nonlinear superposition of a disk 
and a black hole was first obtained by Lemos and Letelier [Jj. Perfect fluid 
disks with halos ^3] and charged perfect fluid disks [20] were also studied. 
For a survey on self gravitating relativistic thin disks, see for instance |21j . 

In the works cited above an inverse style method was used to solve the 
Einstein equations, i. e., the energy-momentum tensor is computed from 
the metric representing the disk. Another approach to generate disks is by 
solving the Einstein equations given a source (energy- momentum tensor). 
This has been used by the Jena group to generate several exact solutions of 
disks HH-HH. 

Even though in a first approximation thin disks can be used as usefull 
models of galaxies, in a more realistic model the thickness of the disk should 
be considered. The addition of a new dimension may change the dynamical 
properties of the disk source, e. g., its stability. Thick static relativistic disks 
in various coordinate systems were presented in [201 • The method used to 
construct the thick disks is a generalization of the well known "displace, cut 
and reflect" method used to generate thin disks from vacuum solutions of 
Einstein equations. This generalization adds a new step and thus can be 
named "displace, cut, fill and reflect" method. In 30 a particular function 
with properties that will be discussed later was used to "fill" the disks. 

In this article we present a class of the functions mentioned above and 
use them together with the Schwarzschild metric to construct more mod- 
els of exact relativistic thick disks. In Sec. [21 we discuss the main idea of 
the "displace, cut, fill and reflect" method and use the Newtonian Kuzmin- 
Toomre disk to put constraints on the parameter in the class of "fill" func- 
tions so that the disks are physically acceptable. Then in Sec. El we take 
the Schwarzschild solution in isotropic cylindrical coordinates to generate 
thick disks. These disks have equal radial and azimuthal pressures but are 
different from vertical pressures. In Sec. the same procedure is repeated 
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(a) (b) (c) 

Figure 1: Construction of a thick disk by the "displace, cut, fill and reflect" 
method. In (a) the space with a singularity is displaced and cut by a plane 
(dashed line). In (b), after disregard the part with singularities, a thick shell 
is put below the plane. In (c), the resultant conflguration is reflected on the 
bottom surface of the shell. 

in Weyl coordinates. The resulting disks have radial tensions that have the 
same modulus as the vertical pressures, azimuthal tensions near the center 
and azimuthal pressures for larger disk radii. In Sec. [3 thick disks are con- 
structed in canonical Schwarzschild coordinates, which were not previously 
studied. Here an additional function must be added in order to generate 
exact disks. The disks show similar characteristics of those constructed in 
Weyl coordinates, but radial tensions are different from vertical pressures. 
Finally, in Sec. El we summarize the main results. 

2 Newtonian Thick Disks 

In order to construct a thick disk we use a modiflcation of the "displace, 
cut and reflect" method, due to Kuzmin [SJ. It can be divided in the 
following steps (Fig. P): first, choose a surface (usually the z = plane) that 
divides the usual space in two parts: one with no singularities or sources, 
and the other with singularities. Second, disregard the part of the space 
with singularities. Third, put a thick shell below the surface. Fourth, use 
the bottom surface of the shell to make an inversion of the nonsingular part 
of the space. The result will be a disk of infinite extension and thickness 
twice of the shell. In the original method the surface that divides the space 
is used to make an inversion, resulting in a space with a singularity that is 
a delta distribution with support on z = 0. Because of the additional step, 
the modification can be named "displace, cut, fill and refiect" method. 
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Mathematically the procedure is equivalent to make the transformation 
z — > h{z) +6, where 6 is a constant and h{^z) an even function of z. In order 
to generate thick disks with well defined properties, the function h(^z) and 
its first derivative should be continuous in the region —a < z < a, where 2a 
is the disk thickness. Furthermore, the first and second derivatives of h{z) 
with respect to z should be chosen such that the mass density of the disk (a) 
is non-negative everywhere and (b) is a monotonously decreasing function 
of r and z. In the case of Newtonian gravity, the potential <I>(r, z) satisfies 
the Laplace equation 

$,rr + — + '^,zz = 0. (1) 

r 

After we make the transformation z h(z) + 6, Eq. leads to 

V^cD = h"<^^h + {h'^ - (2) 

where primes indicate differentiation with respect to z. From Eq. @ we 
have that, if \h'\ = 1 and thus h" = when \z\ > a, the mass density 
vanishes outside the disk. For \z\ < a, we get from Poisson equation 

^ + _ 1)$ . (3) 
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If we start from an even and positive polynomial for h"{z), we arrive at 
a class of functions h{z) that satisfy all the requirements stated above given 
by 

—z + C, z < a, 

h{z) = I Az'^ + Bz'^"-+'^, -a<z<a, 
z + C, z > a, 



(4) 



with 



A 



2n + 1 — ac 



B 



ac 



1 



C 



o(2n + 1 + ac) 



Ana ' 4n(n + l)a2"+i ' " 4(n + 1) 

Here n = 1, 2, . . . , and c is the jump of the second derivative on z = ita. The 
special case when ac = 1 was considered in [HOI- As an example we consider 
the gravitational potential of a mass point in cylindrical coordinates 

Gm 



(5) 



By doing the transformation z — > h{z) + 6 in the previous potential and 
using Eq. Q, we obtain the mass density 



m 



h"(h + b) + h 



12 



i?3 



1 3(1 

- + — 



(6) 
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with i?^ = + (/i + 6)^ and the variables and parameters were rescaled in 
terms of the disk half-thickness: f = r/a, h = h/a, b = b/a, fa = m/a, 
p = c? p and with c = c/a and z = zja'm. Eq. @. The mass density will be 
positive if 

h!'(h + 6) + - 1 > 0. (7) 
When i = 0, condition ^ reduces to 

~2n + 1 - c r 2n 

h 1>0^6> , 8) 

2n - - 2n + l-c ^ ' 

and c < 2n + 1 to make b positive. When z = 0, Q gives c{h + 5) > 0, or 
c > 0. Thus the parameter c is restricted to < c < 2n + 1. 
The total mass M. of the disk is easily calculated 

M. = \ / / pf df dzd(/9 = m, (9) 

Jo Jz=-l Jf=0 

thus the disks have finite mass. 

In Fig. [21 we plot the curves of h{z) and its first and second derivatives 
with n = 1 for c = 0, 0.5, 1 and 2. We also plot the curvature k of h, 
calculated from the expression 

\h"\ , , 

(10) 



(1 + /i'2)3/2 



Fig. Eland Fig.jllshow the mass density for a Newtonian thick disk calculated 
with the function h{z) depicted in Fig. I2ta) and with parameters m = 1 and 
6 = 2. In Fig. Owe plot again h{z) with n = 2 and the same values for c 
as in Fig. |21 and the mass density for a Newtonian disk using this function 
with parameters m = 1 and 6 = 2 is shown in Fig. j^l and Fig. [71 We note 
that the mass distribution along f = has a similar shape to the curves 
for the curvature k. Above some value of the jump c, the maximum of the 
mass density at (f, z) = (0, 0) becomes a local minimum point and two other 
maximum points appear. This may be interpreted as a split from one into 
two disklike distributions of matter, a configuration that does not seem to 
be physically reasonable. Thus we restrict the interval of values for c such 
that the mass density has only a central maximum. The critical points of 
the curvature function are given by 



+ -3/i'/i"2 = 0. (11) 
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Figure 2: (a) The function h{z) and (b) its first and (c) second derivatives 
with n = 1 for c = 0, 0.5, 1 and 2. (d) The curvature k{z) for the same 
parameters. 
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Figure 3: The mass density Eq. © for a Newtonian thick disk with param- 
eters m = l, 5 = 2, n = l and (a) c = 0, (b) c = 0.5. Some levels curves of 
the density are displayed on the right graphs. 
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(b) 




r 



Figure 4: The mass density Eq. for a Newtonian thick disk with param- 
eters m = l, 6 = 2, n = l and (a) c = 1, (b) c = 2. Some levels curves of 
the density are displayed on the right graphs. 
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Figure 5: (a) The function h{z) and its (b) first and (c) second derivatives 
with n = 2 for c = 0, 0.5, 1 and 2. (d) The curvature k{z) for the same 
parameters. 
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(b) 




Figure 6: The mass density Eq. © for a Newtonian thick disk with param- 
eters m = l, 5 = 2, n = 2 and (a) c = 0, (b) c = 0.5. Some levels curves of 
the density are displayed on the right graphs. 
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(b) 




r 



Figure 7: The mass density Eq. for a Newtonian thick disk with param- 
eters m = l, 6 = 2, n = 2 and (a) c = 1, (b) c = 2. Some levels curves of 
the density are displayed on the right graphs. 
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For n = 1, this condition leads to z = and the roots of 

7{c-lfz^+ll{3-c){c-lf2'^+5{3-cf{c-l)z^-8{c-l) + {3-cf = 0. (12) 

To have z = as only critical point, Eq. (|12() should not have real roots. 
A graphical analysis shows that this happens if c ^ 1.46. For n = 2, the 
polynomial equation in question is 

55(c-l)^5^^+65(5-c)(c-l)2z^+13(5-c)2(c-l)i^-320(c-l)52+3(5-c)^ = 0, 

(13) 

which does not have real roots if c ^ 1.81. 

The level curves in Fig. 0] and in Fig. EE also indicate some undesirable 
features of the disks: some isodensity curves show that the density decreases 
towards the disk's center. This happens even for n = 1, c = 1 in Fig. |lja) 
and n = 2, c = in Fig.lHfa). Thus, in practice, the ranges for the parameter 
c that allow physically acceptable disks is even more restricted than those 
stated above. Using graphical analysis, we find that large values for the 
"cut" parameter b allow a larger "good" range for the parameter c, but in 
the next sections when the procedure will be applied in general relativity, 
large values for b will also mean less relativistic disks. 



3 Thick Disks from the Schwarzschild Metric in 
Isotropic Coordinates 

In cylindrical coordinates {t,r, z,(p) the isotropic metric takes the form 



.21/, 



(14) 



where the functions and fi depend on r and z. The Schwarzschild solution 
for metric Eq. (|14j) can be expressed as 



In 
2 In 



2R 



m 



2R + m 
m 



1 + 



2R 



(15a) 
(15b) 



+ z^. The "displace, cut, fill and reflect" method 



where m > and 

applied to the above solution is equivalent to put = r"^ + {h + 6)^ where 
6 > and h[z) given by Eq. @. The continuity of the metric function and 
its derivatives on z = ita is satisfied by the continuity of h[z) and h'{z) 
on z = iba. For \z\ > a, Eq. (|14|) satisfies the vacuum Einstein equations. 
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For \z\ < a, the components of the energy-momentum tensor of the disk are 
computed from the Einstein equations 



Rn 



1 



9abR, 



(16) 



where we use units such that c = SttG = 1. Defining the orthonormal tetrad 

jya^^a^ya^^-"}, with 



1^'^ = 6-^(1,0,0,0), 



'(0,1,0,0), 



^^ = 6-^(0, 0,1,0), 
Z" = —(0,0,0,1) 



(17a) 
(17b) 
(17c) 

(17d) 



and using Eq. (jlGf) . the energy-momentum tensor can be written as 

Tab = (jVaVb + PrWaWb + p.YaYb + p^ZaZb, (18) 

where a = — T/ is the energy density, pr = is the radial stress, which is 
equal to the azimuthal stress p^ = , and pz = is the vertical stress. 
The effective Newtonian density is given hy p = a -\-pr+Pz+Pip- To satisfy 
the strong energy condition we must have /9 > 0, the weak energy condition 
requires u > and the dominant energy condition requires |pr/c| < 1, 
\Pz/cr\ < 1 and \p^/a\ < 1. Using Eq. (|15aj) - (|15b|) . we obtain [32] 



64m 



(m + 2Rf 

128mR 



3{h + 6)^(1 - h") + R'[h" - 1 + h"{h + b)] 



(19a) 



{m + 2R)^{-m + 2R) 
32m2 



Pr = Pifi 



{fh + 2Rf{-m + 2R) 
{m + 2Rf{-m + 2Ry 



where a = a?a, p = a?p, p^ 



3{h + 6)^(1 - h") + R'[h'^ - 1 + h"{h + b)] 



(19b) 

2(/i + 6)2(1 - + R2[h'^ _ 1 + h"(^h + b)] 

(19c) 
(19d) 



a?Pz, and the other dimensionless 



a Pr, Pz 

variables were defined in Sec. EJ 

From Eq. (|19a() - H19c() we can see that if condition ^ is satisfied, we 
have IT > 0, p > and Pr = p<f > (pressures). Thus b > 2n/(2n + 1 — c) 
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Figure 8: (a) Energy density, (b) effective Newtonian density for a thick 
disk in isotropic coordinates. Parameters: m = l, b = 2, n = l and c = 0. 
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(a) 




Figure 9: (a) Radial and azimuthal pressures, (b) vertical pressure for a 
thick disk in isotropic coordinates. Parameters: m = l,6 = 2,n=l and 
c = 0. 



15 



and < c < 2n + 1. We also have vertical pressure. To ensure nonsingular 
behaviour of the expressions, we impose that < m < 26. 

In Fig. ISl and Fig. |^ we plot the surfaces and level curves for the energy 
density, effective Newtonian density, radial and azimuthal pressures and 
vertical pressure for a thick disk with parameters m = l, b = 2, n = l and 
c = 0. The ratio between any pressure and energy density is less than 0.15. 
Thus all energy conditions are satisfied. The shape of the level curves also 
show that for these parameters the disk is physically acceptable. 



4 Thick Disks from the Schwarzschild Metric in 
Weyl Coordinates 

The general metric for a static axially symmetric spacetime in Weyl's canon- 
ical coordinates {t,r, z,ip) is given by 



(20) 



where $ and A are functions of r and z. The Einstein vacuum equations for 
this metric reduce to the Weyl equations ^21 0^ 



^-.rr H + '^,zz 

r 

A, 



0, 



In these coordinates, Schwarzschild solution assumes the form [,VA\ 

^Ri + R2-2m' 



A 



Ri + R2 + 2m 
{Ri + Rjf - 4m^ 

4:RlR2 



(21a) 

(21b) 
(21c) 



(22a) 
(22b) 



where Rf = r'^ + (m + z)^, R^ = r"^ + {—m + z)'^ and m > 0. By applying the 
transformation z h{z) + 6 on Eq. iP^ - l|22b|) . using Eq. ipTa|) - ipTc|) and 
Einstein equations Eq. (|16() . the energy-momentum tensor can be written as 



Tab = eVaVb + PrWaWb + PzYaYb + P^ZaZ^, 



(23) 
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whith the orthonormal tetrad {V"", VF",y", Z"-} defined as 



F'^ = e-*(1, 0,0,0), 



= e*-^(0, 0,1,0), 



-(0,0,0,1). 



(24a) 
(24b) 
(24c) 

(24d) 



In Eq. (|23|1 e = — T/ is the energy density and Pr = T^, Pz = = —T^, 
Pip = Tip are respectively, the radial, vertical and azimuthal stresses. The 
effective Newtonian density is given hy p = e + pip. The explicit expressions 
are 



g2(*-A) 



|/t"(i?i - R2)RIrI + (1 - h'^) {h + h){Rl - Rl) - m{Rl + rI 



1 



2mP(i?i +i?2) 



RiR2[{Ri + R2Y -^m^] 
{Rl - R2)RlRl - 2r^{Rl - R^) 



+ 



e^^^-^\l-h'^){Ri- R2) 

m\R\ 



2mpf\Ri + R2) 



+ 



Pz 



RiR2[{Ri + R2? - ^m^] 

{Rl - R2)RiRl - 2r^{Rl - R^) 

4m\n^-^){h'2 -l){h + lf 

i?^i?|(i?l + i?2)^ 
_ 47fi2e2('f-A)(l_/i'2)(/j + 5)2 

i?^i?2(-Rl + -^2)^ 



(25a) 

(25b) 

(25c) 
(25d) 

(25e) 



From Eq. (|25d|) and Eq. H25e() we have vertical pressures and radial tensions. 
From Eq. ()25b|) the condition e > on |z| = 1 gives 



2rhr'^{Ri + R2 



RiR2[{Ri + R2y - irh^] 
Since Ri > r, R2 > f, and Ri + R2 > 2b, we have 

mb 



< 1. 



(26) 



fj2 _ ^2 



< 1, 



(27) 
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(a) 




r 



Figure 10: (a) Effective Newtonian density, (b) energy density for a thick 
disk in Weyl coordinates. Parameters: m = l, 6 = 2, n = l and c = 0. 

which is equivalent to m < (-v/S — l)b/2 It is not easy to obtain 

other constraints over the parameters m and b in order to satisfy ah energy 
conditions. The analysis is better done graphically. 

In Fig. ^1 and Fig. ^2 we plot the surfaces and level curves for the effec- 
tive Newtonian density, energy density, vertical pressure and the azimuthal 
stresses for a thick disk with parameters rh=l,b = 2,n=l, and c = 0. For 
these values we have p > and e > everywhere. The ratio between any 
pressure or tension and energy density is less than 0.35, thus the dominant 
energy condition is also satisfied. 
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(a) 




Figure 11: (a) Vertical pressures, (b) azimuthal stresses for a thick disk in 
Weyl coordinates. Parameters: m = l, b = 2, n = l and c = 0. 
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5 Thick Disks from the Schwarzschild Metric in 
Schwarzschild Coordinates 



The Schwarzschild metric in Schwarzschild canonical coordinates (t, r, 9, ip) 
is written as 



ds' = -{l-^] dt' + T^^^dr^ + r\de' + sin2 Od^'), (28) 



r 



where m is a positive constant. In cylindrical coordinates {t, R, z,<p), Eq. 
H28|) can be cast as 

2 / 2m\ 2 Tl^ - 2mz'^ 2 4mi?z 



ds = - 1 - — dt^ + — — -dR^ + -TT-- -di?d2+ 

V n J 7^2(7^ - 2m) n^{n - 2m) 

where 7^^ = R? + . The "displace, cut, fill and refiect" method cannot be 
directly applied to metric Eq. 1)29^ . since the component is not an even 
function of z. But we can multiply it by an odd function of z, which we 
choose to be the first derivative of h[z). So we make the transformations 
z h{z)+b and qrz — > 9Rzh'{z) in Eq. ((^ . and use the Einstein equations 
Eq. H16|) to compute the components of the disk's energy-momentum ten- 
sor. With h{z) defined by Eq. @, the metric Eq. H29|) satisfies the vacuum 
Einstein equations for |z| > a and the metric functions together with their 
derivatives with respect to z are continuous on 2; = ita. Note that the func- 
tion multiplying g^i^ niay be other kind of function not necerrarily related 
to h[z)^ but its possible form is also limited by the requirements that the 
generated disks are physically acceptable, as was discussed in Sec. I2j The 
physical variables of the disk are obtained by solving the eigenvalue problem 
for 

Tl^i^ = XC- (30) 
We find that T"'^ can be put in the form 

ri^ab ^ ^jjajjb ^ p^yayb ^ p X'^X'' + p^WW^ (31) 
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where 



-T*, = ^^(1,0,0,0), 

V-9tt 



p± = ± ^V(7T-r|)2 + 4r;r/, 

= (o,y,y^o), X" = (o,x^x^o), 

p^ = r^, M^« = 1(0,0,0,1), (32) 



and 



1 T'' - B 

X'^ = , = r = / ~ ■ (33) 

The effective Newtonian density reads p = e+p++p_+P(^ = e+T^+T^+p^. 

Although the expressions are exact, we do not state them explicitly, since 
they have dozens of terms. The analysis must also be done graphically. We 
find that is a pressure along a direction mostly vertical and p- is a 
tension along a direction mostly radial. In Fig. 11211111 we graph the surfaces 
and level curves of the energy density, the effective Newtonian density, the 
pressure mostly vertical, the tension mostly radial and the azimuthal stresses 
with parameters m = l, 6=5, n = l and c = 0. The variables and 
parameters were rescaled as in the previous sections. While densities and 
the radial tension decrease monotonously with f and z — > ±1, the upward 
pressure first increases with radius and then decreases monotonously. The 
azimuthal stresses are negative near the disk's center, then increase to a 
positive maximum at the z = plane and then decrease monotonously. In 
some sense this new class of disks has similar characteristics of those obtained 
in Weyl coordinates, but here the almost radial tensions are different from 
the almost vertical pressures. Unfortunately the dominant energy condition 
does not hold with respect to the azimuthal stresses. Fig. I14f b) shows the 
level curve of \pip/e\ = 1. For r ^ 8 the dominant energy condition is not 
satisfied. This also seems to happen for larger values of b. In Fig. I14r c) 
we plot the lines of flow calculated by numerically solving the differential 
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(b) 




r 



Figure 12: (a) Energy density, (b) effective Newtonian density for a thick 
disk in Schwarzschild coordinates. Parameters: m = l,6 = 5, n = l and 
c = 0. 
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Figure 13: (a) The pressure mostly vertical, (b) tension mostly radial for 
a thick disk in Schwarzschild coordinates. Parameters: m = 1, 6 = 5 and 
c = 0. 
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Figure 14: (a) The azimuthal stresses, (b) the level curve of \p^p/f\ = 1 and 
(c) the lines of flow calculated from Eq. (|34ap - H34b|) . Parameters: m = 1, 
6 = 5 and c = 0. 
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equations 



dX-" X"- 



(34a) 
(34b) 



The vertical and horizontal lines in Fig. I14f c^ are, respectively, the lines 
of flow associated to Eq. (|34a|) and Eq. H34b|) . We note that the lines are 
very parallel to the coordinate axes. This happens because the component 
is about two orders of magnitude smaller than the other components 
of the energy-momentum tensor, thus the eigenvalues p± are almost equal 
to the corresponding diagonal elements. We call then the eigenvalue p+ 
with eigenvector V"" the pressure mostly vertical and the eigenvalue p- with 
eigenvector X"" the tension mostly radial. 



6 Discussion 



The "displace, cut, fill and reflect" method and a class of functions used 
to "fill" was used to construct new classes of thick general relativistic disks 
that generalize the models studied in a previous work jSU]. The Newtonian 
Kuzmin-Toomre disks were used to put constraints on the parameter in the 
class of "fill" functions so that the disks were physically acceptable. Then 
the Schwarzschild solution was used to construct disks in isotropic cylindri- 
cal coordinates, Weyl coordinates and canonical Schwarzschild coordinates. 
In isotropic coordinates the disks have equal radial and azimuthal pressures 
(isotropic fluid) but different vertical pressures. Disks in Weyl coordinates 
present radial tensions that are equal in modulus as vertical pressures, az- 
imuthal tensions near the disk's center and azimuthal pressures for larger 
radii. All disks are in agreement with the strong, weak and dominant energy 
conditions. 

In canonical Schwarzschild coordinates the "displace, cut, fill and reflect" 
method cannot be applied directly and an additional odd function must be 
added to generate exact disks. We find that the disks have tensions along a 
direction mostly radial that are different from the pressures along a direction 
mostly vertical, and the azimuthal stresses have similar behavior as in disks 
in Weyl coordinates. These disks are in agreement with the strong and weak 
energy conditions, but the azimuthal stresses do not satisfy the dominant 
energy condition. 
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